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ABSTRACT 

I demonstrate that an effect similar to the Romer delay, familiar from timing radio pulsars, should 
be detectable in the first eclipsing double white dwarf (WD) binary, NLTT 11748. By measuring the 
difference of the time between the secondary and primary eclipses from one- half period (4.6 s), one 
can determine the physical size of the orbit and hence constrain the masses of the individual WDs. 
A measurement with uncertainty < 0.1s — possible with modern large telescopes — will determine the 
individual masses to ±O.O2M0 when combined with good-quality (< lkms~^) radial velocity data, 
although the eccentricity must also be known to high accuracy (±10"'^). Mass constraints improve as 
(where P is the orbital period), so this works best in wide binaries and should be detectable 
even for non-degenerate stars, but such constraints require the mass ratio to differ from one and 
undistorted orbits. 

Subject headings: binaries: eclipsing — stars: individual (NLTT 11748) — techniques: photometric — 
white dwarfs 



1. INTRODUCTION 

Sin ce the discovery of binary pulsars (jHulse fc Tavloii 
|1975[ ). precision timing (typical uncertainties < 1 /is) 
has been used to deriv e a variety of physical con straints 
(see the discussion in iLorimer fc Kranieil 120041 ) . The 
arrival-time delay across the orbit (the Romer delay ^) 
immedi ately gives the projected sem imajor axis of the 
pulsar (Bland ford &: Teukolskvl 119761 ). This, especially 
when coupled with the relatively narrow mass distr i- 
bution of neutron stars (|Thorsett fc Chakrabartviri999l ). 
constrains the mass of the companion. 

I contrast this with eclipse timing of planetary sys- 
tems (typically uncertainties are > seconds). Here, with 
a mass ratio ~ 10^'^ the radial velocity curve gives 
a limit on the mass of the companion planet. With 
transiting systems, smi fs 1 and the mass of the 
planet is further constrain ed but not known uniquely 
(jCharbonneau et all l2000t ). although with knowledge 
of the stellar paramet ers one can infer t he planetary 
mass and radius (e.g., iBrown et al.ll200i|) . If individ- 
ual eclipses can be timed to high precision (and here I 
mean both primary and secondary eclipses, i.e., transits 
and o ccultations) , one can learn more about the system 
(e.g., I Wind [201 Of ). Variations in the eclipse times can 
unvei l the presence o f additional bodies in th e system 
(e.g.. lAgol^et al l 120051: THolm an fc Murravll2005D . preces- 
sion (e.g^Miralda-Escudc 20021 ). and kinematics of the 
system (|Rafikov.i2009 ) . 

With the recent discovery (jSteinfadt et all [201Qb( ) of 
NLTT 11748, an eclipsing double white dwarf (WD) 
binary with a tight enough orbit that the binary will 
merge within a Hubble time, a whole new series of ques- 
tions may be asked. The initial constraints are the ra- 
dial velocity amplitude of the lighter object (owing to 
the inverted mass-radius relation of WDs, this object 

1 Hubble Fellow 

^ The Romer delay is named after O. Romer, who used deviations 
from periodicity in the eclipses of lo by Jupiter to deduce a finite 
speed of light jSterken 2005b. ). 



is the larger and brighter member of the system) and 
the widths and depths of both transit and occultation. 
From t his, assuming a cold C /0 WD for the heavier 
object, iSteinfadt et al.l ()2010b[ ) were able to limit the 
masses and radii of both objects, but could not determine 
unique constraints. Measurement of spectral lines from 
the fainter object would determine both masses uniquely, 
but this is challenging as the fainter object is only 3.5% 
of the flux of the brighter. 

A number of other close WD binaries have been 
discovered in the last 2 years (see Table [T] for those 
with undetermined inclinations). Most of them, like 
NLTT 11748, appear to have a low-mass (< 0.2 M©) 
He WD in orbit with a more massive (0.5-1.0 Mq) C/0 
WD. Such systems are of interest because of their even- 
tual evolution, with mass transfer brough t on by grav- 
itational radiation (Nclcma ns et al.l I2001D and are pre- 
sumed to be the progenitors of highly variable objects: 
R CrB stars, AM CVn binaries, and Ty pe la supernovae 
(|Iben fc Tutiikovl[T98l I WebbinkI IToM l . Many of these 
binaries are also of immediate interest as verification tar- 
gets for the Laser Interferometer Space Antenna (LISA) 
mission (Nclcmans 2009). 

Of the 11 compact WD binaries known, only 
NLTT 11748 is known to be eclipsing, but searches for 
the other sources are not uniformly constraining and ad- 
ditional systems may yet be discovered. The flux ratios 
vary for the systems, and in some cases it may be easier to 
directly measure the radial velocity curves for both mem- 
bers of the binary. Without two radial velocity curves, 
mass constraints are limited. Such constraints are invalu- 
able in understanding the detailed formation histories 
and expected evolution of these systems as well as in de- 
termining the mass-radius relation from eclipse measure- 
ments. Moreover, their use as LISA verification sources 
is improved by accurate knowledge of the binary parame- 
ters. In this Letter, I discuss an effect that uses precision 
timing of the eclipses in such double WD systems to help 
constrain the individual masses of the WDs. This tech- 
nique is known in other contexts, being common in ra- 
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TABLE 1 

Double WDs That Will Merge Within a Hubble Time And May Be Eclipsing 



Object 


^brb 


K2 


Ml 


Ri 


Ma 


R2 




Refs. 




(hr) 


(kms-i) 




(Re) 


(Mq) 


(Re) 


(s) 




SDSS J1053+5200 


1.02 


265 


0.26 


> 0.017 


0.20 


0.04 


0.2 


1,2 


SDSS J1436+5010 


1.10 


347 


0.45 


0.014 


0.22 


0.04 


0.7 


1,2 


SDSS J0849+0445 


1.89 


367 


0.65 


0.012 


0.17 


0.05 


2.0 


1 


WD 2331+290 


4.08 


156 


0.39 


0.015 


0.32 


0.016 


0.5 


3,4,5 


SDSS J1257+5428 


4.55 


323 


0.92 


0.009 


0.15 


0.04 


4.7 


6,7,8 


NLTT 11748 


5.64 


271 


0.74 


0.010 


0.15 


0.04 


4.6 


9,10 


SDSS J0822+2753 


5.85 


271 


0.71 


0.010 


0.17 


0.04 


4.7 


1 



Ref erences. — 1: 'Kilic ct alj ||20TOD : 2: [M ullally ct al.l 

l[200gh : 3: [Marsh. Dhillon. fc Duc k (199^ 4 : [NeTemans et a l. (2003); 

5; 'Liebert. Bergeron, fc Holberd li2005i) 6j IBadenes efall ll2009r'): 7: 

[KuTk arni & van Kerkwiik (2010); 8: IMarsh et al.l 120101 1: 9: IKawka fc VennesI II2009I '): 10: 
l^temfadt et aLi t2010b ') 

Note. — The values for Mi assume an edge on orbit, i.e., sin i = 1. The values for 
Rl w ere calculated from those values assuming a cold C/O WD lAlthaus fc Benvenutol 
1199811 . Additional double- WD binaries exist, but they have inclinations known to exclude 
eclipses. 



dio pulsar systems a nd planetar y systems ( [ Knutson et aLl 
[S)07; Hcbr ard et al . 2010; A gol et al.l[2010( ). although in 
the latter it is largely a nuisance parameter and does 
not constrain the systems. I discuss its applicability to 
eclipsing double WD systems, the required observational 
precision and the resulting accuracy. 

2. LIGHT TRAVEL DELAY AND MASS CONSTRAINTS 

In a system with a circular orbit, one often speaks of 
the primary and secondary eclipses as occurring exactly 
1/2 period apart, but this is not the case. If the members 
of the binary are of unequal mass the finite speed of light 
will cause an apparent shift in the phase of the secondary 
eclipse from P/2, where P is the period of the binary 
(|Loebll2005l: iFabrvckvllMol ) . This is similar to the shifts 
in eclipse timing caused by a perturb ing third body on 
a binary system ([Schneider & Dovl(^ | 1995l: iDovle et alJ 
Togs'; 'Dee g et al l \200(k iSterkeni i2005al: iLee et al.l 120091: 
Qian ct al.l l2009f ). although here one only requires two 
bodies and the frequency of the shift is known. 

In the case of a planet with mass m ^ M orbiting a 
star with mass M, one has a primary eclipse when the 
planet is in front of the star. The light is blocked at time 
t = Q. However, that light was emitted earlier by the 
star, at time ti — — a/c, since it traveled a distance 
a (the semimajor axis). For the secondary eclipse, the 
light is emitted by the planet at time t = P/2 but is 
blocked a/c later, at t2 = P/2 + a/c. The difference of 
these times exceeds P/2 by At^T — t2~ti — P/2 = 2a/ c, 
the sought-after quantity. 

For two finite masses, I consider two objects orbiting 
their center of mass with period P, masses Mi and M2, 
and semimajor axis a. The total mass of the system is 
M ^ M1+M2, and of course 4n^a^ = P^GM; the first 
object orbits at a radius ai — a{M2/M) and the second 
object orbits at a radius 02 = a{Mi/M). 

Near primary eclipse, the primary is at [x,y\ = 
[2T:ait/P, ai] and the secondary is at [—2'Ka2t/P, —0.2] at 
time with the observer at [0, — cx)]. I project the image 
of the two objects to the barycenter at y = 0. This gives 
xb.i — 2'Kai{t ~ ai / c) / P and 2:3.2 = — 27ra2(t + a2/c)/P. 
Eclipses occur when these are equal, which has the so- 
lution ti = (ai — a2)/c. Near secondary eclipse, the 



primary is at [—2nai{t — P /2) / P,ai\ and the secondary 
is at [2iTa2{t — P/2)/P, — 02]. Following the same argu- 
ment, eclipses occur when t2 — P/2 = (02 — a\)/c. So 
the eclipses differ by ^2 — = P/2 -I- 2(a2 — ai)/c. The 
light-travel delay is again A^lt = ^2 — ^1 — P /'2^ 

reaching 2a /c when M2 <C A/i, as expected. 
From Kepler's laws the mass function K\P /2'nG = 
Aff sin^ iM~^, where K2 is the radial velocity ampli- 
tude of object 2, and since it is a transiting system, 
sin I K, 1. Substituting for a and the masses, the time 
delay in terms of observables and the mass ratio q (where 
q = M2/M1 < 1) is: 

PK9 , 

AtLT = ^1-9- 2 

TTC 

2.1. Magnitude and Detectability 

The eclipse duration for a cir cular orbit is roughly 
T « 2i?2-P/(27ra) « 3 minutes (|Wh^ IMol) : the du- 
ration of ingress/egress r is decreased by a factor of 
R1/2R2 « 8, T « RiP/{2TTa) fa 20 s (numerical results 
are for NLTT 11748); and ingress/egress are sharpest at 
inclinations of exactly 90°. The accuracy of the eclipse 
time determination largely depends on the duration of 
the ingress/egress and the total number of photons ac- 
cumulated during ingress/egress, since the bottom of the 
eclipse is only slightly curved (for the primary eclipse) 
if not flat (second ary eclipse), and one can derive (e.g., 
ICarter et al.l[2008h 

where the eclipse has fractional depth 6, each observation 
has fractional uncertainty a, and there are A'obs obser- 
vations during r. This holds in the limit that the noise 
is un correlated (cf. ICarter fc Wiiml 120091: iSvbilski et all 
l2010f ). which should be true at the level discussed here 
(photometric precision of > mmag). So, at^ scales as the 
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duration of ingress/egress divided by the total signal-to- 
noise ratio accumulated during that portion of the or- 
bit. Since A'obs oc r for a constant observing cadence, 
at^ oc 77. A star with V = 16.5 mag like NLTT 11748 
gives roughly 0.1 photon s~^cm~^ or 2 x 10^ photons de- 
tected during a 20 s ingress/egress with a 4 m telescope. 
For an eclipse depth of 5% this means a precision on 
individual eclipse times of < Is. 

For general binary systems, I can rewrite the expres- 
sion for AtLT in terms of the primary mass Mi, q, and 
P (eliminating K2): 



AtLT 



/2GA/iP2 



1/3 



(i-g) 

(l + g)2/3 



(4) 



For systems with primaries that are typical C/0 WDs 
and with secondaries that are He WDs, with Mi = 0.5 — 
IMq and q — 1/6 — 1/2, AtLT goes from 0.5 s to 7 s for 
periods of 0.5-lOhr (Tabled]). 

Eclipse depths are functions of the radii and tempera- 
tures of the WDs, as well as the bandpass, and are hard 
to predict with any generality. The main factor that will 
change systematically for other double WD systems is r 
itself, which is oc RiP^^^. This means the signal-to-noise 
ratio (i.e., detectability) is /S.tij^/(Jt^ oc \J PjRi, so the 
effect is easiest to see in long-period binaries. As the 
mass ratio approaches 1 the magnitude of the delay de- 
creases, limiting its utility, but in such systems it may 
be easier to search for the second set of spectral lines 
(depending on the temperatures of the objects). 

This effect should also be present in partially degener- 
ate (sdB+WD) or non-degenerate binary systems. The 
precision on the eclipse times goes as y'r oc V^smaii- 
If the binary is wide enough that the larger star(s) 
are undistorted by tides and hence the orbit remains 
strictly periodic, the overall detectability AtLT/ct^ oc 
yj Pj i?smaii Can actually increase over the double WD 
case I have been considering. Requiring a oc i?iarge to 
minimize tidal distortions, which scale as (i?iarge/a)'^, 
for a system with a O.2M0 M star one needs periods 
of > 1 day to have tidal effects that are as small as in 
the double WD systems. For the ingress/egress duration 
the radius of the smaller object i?smaii increases from 
^ O.Oli?© to ^ 0.1i?o (for a O.IMq M star companion, 
for example), so if the period increases by more than 
a factor of 10 then the wide system is more easily de- 
tectable (the probability of eclipse does decrease as i?/a, 
though, and both primary and secondary eclipses must 
be seen). However, the ephemeris must be known suffi- 
ciently well with tight enough limits on (or measurements 
of) eccentricity so that the light-travel delay is the only 
deviation from regularity (see below). 

For NLTT 11748, I recognize that K2 is the ra- 
dial velocity that was measured since the heavier ob- 
ject is the fainter one. So, q « 0.15/0.71 = 0.21, 
K2 = 271 kms~^, and P = 5 .64hr, which give 
AtLT = 4.6 s. iSteinfadt et al.l (|2010b( ) measured individ- 
ual eclipse times to ^ 10 s, making it hard to detect an 
effect like this. However, this was using 45 s exposures 
on a 2 m telescope, while the ingress/egress duration 
was only 20 s. Increasing to 4 m or 8 m will improve 
the S/N of individual exposures by a factor of 4-16, and 
using a cadence better matched to the orbit will help 



as well, driving eclipse time uncertainties to < 1 s (as 
above). This is sufficient to detect A^lt; below I discuss 
how well one can measure it and what constraints one 
can get from it. 

2.2. Comparison With Eccentricity 

The above discussion considered circular orbits. For 
eccentricity e > the situation changes. I note that 
the objects in Table [T] have orbits that are consistent 
with circular orbits, although quantitative limits for e 
are not always given. This follows from their expected 
evolutionary histories, where common-envelope evolution 
(|Nelemans et al.l 120001) should have circularized orbits. 
Nonetheless, in case our understanding of these systems 
is incorrect or some further evolution (such as interaction 
with another body) may have caused non-zero eccentric- 
ity, I consider the effect of a non-zero eccentricity on our 
detection of AtLT- 

First, there a re changes to the expression for AtLT 
(|Fabrvckvll20ll : 



At 



LT 



(AtLT)e=0X 



1 



1 — sin uj 



(AtLT)e=ox(l-e2 



(5) 

where u is the argument of pericenter. As (AtLT)e=o 
is small to begin with, this is unlikely to be signifi- 
cant. More important, though, is that an additional term 
changes the relativ e timing of th e primary and secondary 
eclipses. Following IWinnI ()2010[ ) : 



Atp 



Pe 



■ cos w 



(6) 



and the primary eclipse also changes duration relative to 
the secondary eclipse by the ratio 1 + e sinw. To compare 
AtLT and Ate means effectively comparing i^2(l — 9)/c 
and ecoscj. For if 2 ^ 300kms~^, this means that one is 
sensitive to e ~ 10"'^ (although a; is poorly determined 
for low e). If it can be asserted for some independent 
reason (i.e., evolutionary assumptions) that e ^ 10"'^ 
then one can treat any measured At as coming from light- 
travel delay. But if not one must be more careful. 

Fortunately, eccentricity can be constrained from the 
radial veloc i ties. Adopting the small-e limit as in 
ILange et al.l (|2001D . spectra can limit the eccen- 

tricity to CTe ~ 2av / K2y/ N-RY , where is the preci- 
sion of the individua l velocity measurements (see also 
iGaudi fc Whmll2007D . With > 100 observations with 
< lkms~^ precision the eccentricity can be limited (in- 
dependent of a;) to <C 10"'^, and hence can identify 
whether any measured time delay has a contribution from 
an eccentric orbit. This requires dedicated radial velocity 
measurements over one or more full orbits, but is achiev- 
able with current instrumentation. At this level one must 
also account for additional effects such as light-travel 
delay in the spectroscopic analysis (jZucker fc Alexandeii 
l2007f ). Tidal disto rtions can a lso mimic eccentricity in 
radial velocity fits (|Eatonll2008l ). but these can be identi- 
fied photometrically and are expected to be quite small, 
^ lO^'*, except in the most compact systems. 

2.3. Mass Constraints 

Equation (j4]) gives an independent constraint on the 
mass ratio q, which helps break the degeneracy in the 
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Fig. 1. — Constraints on individual masses from a single radial 
velocity constraint K2 and a light-travel delay AtLT- I show the 
constraints from each measurement individually (along with 3o" 
ranges) as the diagonal lines. The contours show Icr, 2(t, and 3(j 
joint confidence contours on Mi and A/2; their covariance is appar- 
ent. This system has parameters similar to those of NLTT 11748, 
and I assumed ax = Ikms"'^ and (ja = 50 ms, which is rather 
optimistic. 

mass function to measure the masses of the stars indi- 
vidually. For the individual masses 



Ml 



M2 



2ttGP 



K2 



AiLTTTc) 



2nGP 



^] (7) 



Assume that I measure K2 ± <Jk and A^lt ^ <^a {P 
is typically known to much higher precision); I also as- 
sumed e = 0. How well can I determine the individual 
masses? I know q to: 



{Rial + Atlr^al) 



(8) 



I now wish to see with what precision I can estimate 
the masses from the observations. Doing standard error 
propagation, 



dMi _ fPA'f 
dMi 



V 27rG 



1/3 



q 



dAt 



LT 



\ P^G 



1/3 



■q) 



(1 + 9) 



(9) 



These are the contributions of the uk and cta to the 
uncertainty on the mass, i.e., a\i — aj^ldM/dKl^ + 
(j\ \dM / dAtin:\ ■ For the two terms to be compara- 
ble requires cta — ((<7 + h)P/2i:c)<JK ~ lOPhrC/s-.kms nis 
(where Phr is the period in hr and aK^kms is the un- 
certainty on if 2 in kms~^). With those, I would have 
ctmi ~ O.OIM0 for periods P > Ihr and mass ratios 
q ^ 0.25. The constraint on M2 is similar. I illus- 
trate this in Figure [TJ where I show mass constraints 
on NLTT 11748 from hypothetical time-delay measure- 
ments. 

However, it is likely that the uncertainty from the ra- 
dial velocity amplitude will be considerably less than that 
from Atif£: individual velocity measurements can eas- 
ily have uncertainties of a few kms~^ with a large tele- 



scope, and combining enough of them to give a mean- 
ingful constraint on the eccentricity will likewise end up 
with gk = ctv/VNbv < lkms~^. Getting a compara- 
ble constraint on the time delay seems implausible: with 
individual times measured to ~ Is, > 10* eclipses are 
needed to get a a < 10 ms, and only one time delay is 
measurable per orbit. This means one needs significantly 
higher signal-to-noise per observation than the several 
hundred I have been assuming here or that I will be lim- 
ited by a A- In this case, the precision on Mi improves 
as a Ml ~ P~'^^^aA] including the difficulty in detecting 
the delay a a ^ P^^^, the determination of Mi improves 
as P~^/'^. The only ways that long periods are penalized 
are in terms of observing strategy, as for long periods the 
time to get enough eclipses measured will grow long as 
well, and for the probability of detecting an eclipse in the 
first place since that decreases as 1/a. 

As shown in Figure [U the joint probability distribution 
for (Mi,M2) is strongly correlated between Mi and M2, 
with a much stronger constraint on M2 — Mi than on 
M2 + Ml. However, this would even be true — although 
to a lesser degree — if the other radial velocity amplitude 
Ki were measured, especially if it were at lower preci- 
sion because it is much fainter. It is straightforward al- 
though tedious to compute the linear combinations of 
the masses that minimize/maximize the variance which 
would be preferred for fitting. While not perfect, with 
this constraint one would have a much better picture of 
the system. 

3. CONCLUSIONS 

Motivated by the recent discovery of NLTT 11748, the 
first eclipsing double WD binary, I have examined a phe- 
nomenon that affects precision eclipse timing of such a 
system. With knowledge of the individual masses, one 
can put much stronger constraints on the radii of the 
two WDs, the evolutionary history of the system, and 
its expected outcome, not to mention WD atmosphere 
models and models for the interiors of He WDs (e.g., 
iPanei et all [20071 : iSteinfadt et al.ll2010al ). This effect, a 
delay between perfect phasing of primary and secondary 
eclipses, can be used to constrain the individual masses 
of the binary, something difficult to do otherwise. 

I find that the light-travel delay should be detectable 
in the case of NLTT 11748, and possibly in some other 
similar binary systems should they prove to have both 
primary and secondary eclipses: long periods are favored 
both for detecting AiLT and for using it to constrain 
the masses, although long periods do not favor detecting 
eclipses to begin with. The constraints on the individual 
masses can approach ztO.OlAfQ for plausible data-sets, 
and will likely be limited by the precision of the eclipse 
timing, suggesting that an intensive timing effort on large 
telescopes is worthwhile. Detection of a second radial ve- 
locity amplitude would over constrain the system, lead- 
ing to even tighter determinations of the masses. In non- 
degenerate eclipsing binaries, such as those that Kepler 
may discover, the delay should also be detectable for sys- 
tems with orbital periods of greater than a few days, al- 
though it requires that the mass ratio differs from one 
and that no other unmodeled orbit variations be present 
to a high degree of confidence. 

The orbits of the binary members can also be per- 
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turbed by other bodies in the systems, either on shorter 
(planets or other small bodies in close orbits) or longer 
timescales (a distant body) . In both cases, perturbati ons 
in transit timing may be visible (see Ag ol et 311120051 for 
a detailed discussion) . Given the wide variety in possible 
situations it is out of the scope of this Letter to con- 
sider, but any perceived variation in transit timing must 
be compared against the possible presence of additional 
bodies. There could also be effects that alter the per- 
ceived primary versus secondary eclipse times without 
altering the orbit, su ch as hot s pots due to irradiation 
(|Knutson et alJl2007t lAgol et al.ir201Gh or accretion. For 
the former, I note that the incoming radiation in the dou- 
ble WD systems is typically very small, ~ 10~'^ of the 
outgoing radiation. As for accretion, both WDs are well 
inside their Roche lobes and so none is expected. 



I find the fortunate coincidence that NLTT 11748, the 
one object known to be eclipsing, also has the binary 
parameters that lead to the highest value of A^lt among 
similar double WD binaries. Hopefully, with dedicated 
observing this effect will be detected and can constrain 
the NLTT 11748 system even more than is possible today. 
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